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ABSTRACT 

We analyse the possibility that our Universe could be described by the model recently proposed by 
Melia & Shevchuk (2012), where the Hubble scale Rh—c/H is at all times equal to the distance 
ct that light has travelled since the Big Bang. In such a model, the scale factor is proportional to 
cosmic time and there is neither acceleration nor deceleration of the expansion. We first point out 
problems with the very foundations of the model and its consequences for the evolution of the 
Universe. Next, we compare predictions of the model with observational data. As probes of the 
expansion we use distance data of supernovae type la, as well as Hubble rate data obtained from 
cosmic chronometers and radial baryon acoustic oscillations. We analyse the redshift evolution 
of the Hubble parameter and its redshift derivatives, together with the so-called O m diagnostic 
and the deceleration parameter. To reliably estimate smooth functions and their derivatives from 
discrete data, we use the recently developed Gaussian Processes in Python package (GaPP). Our 
general conclusion is that the discussed model is strongly disfavoured by observations, especially 
at low redshifts (z < 0.5). In particular, it predicts specific constant values for the deceleration 
parameter and for redshift derivatives of the Hubble parameter, which is ruled out by the data. 
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1 INTRODUCTION 

The fact that the Universe is currently undergoing a phase of ac- 
celerated expansion has been known for more than a decade. The 
first evide nce came from obs erva tions of distant superno vae type la 
(SNela) bv lRiess et ail Jl998t) and lPerlmutter et al.l dl999h . This result 
has been confirmed by various observations, e.g. of the cosmic mi- 
crowa ve background and its angular power spec trum dKomatsu et alj 
1201 ll) . bar yon acoustic oscilla tions (BAO, e.g. Percival et al. 20ld) 
and SNela JSuzuki et ap2012h : see lWeinberg et all d2012l) for a gen- 
eral review. 

However, the picture that emerges based on these and multiple 
other observations is far from being theoretically understood. In a ho- 
mogeneous and isotropic Universe and assuming general relativity to 
be correct on cosmological scales, one needs some form of dark en- 
ergy as the driving force of the observed acceleration. The simplest 
model that is consistent with observations is the Lambda-Cold Dark 
Matter (ACDM) - a model where the dark energy has a constant en- 
ergy density. 

This theory has no physical understanding at the moment and in 
recent years considerable effort has been devoted to develop and test 
alternative explanations of the accelerated cosmic expansion. These 
efforts include dark energy models with an equation of state that is 
more complex than that of ACDM, inhomogeneous models and the- 
ories assuming deviations from general relativity, known under the 
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collective name of modified gravity [see Copeland et al.l ( |2006|) for a 
review]. 



2 THE R h = ct MODEL AND ITS PROBLEMS 



In a recent paper, iMelia & Shevchukl 1 120 12l hereafter MS 1 2) have of- 
fered a qualitatively different cosmological model, which was meant 
to solve some of the problems of the ACDM. Namely, they used the 
general relativistic framework of Friedman-Lemaitre models with the 
Robertson- Walker metric and, having applied the Weyl postulate, pro- 
posed a model in which there is neither acceleration nor deceleration 
either now or in the past. It is claimed that this is the only viable 
cosmological model as all other models allegedly violate the Weyl 
postulate. In this eternally coasting model, the scale factor grows lin- 
early with time. Additionally, the MS 12 model predicts the Hubble 
parameter to be always exactly equal to the inverse of the age of the 
Universe, i.e. H(z) — t~ 1 (z), where t is the cosmic time and z is 
the cosmological redshift. MS 12 argue that the best-fit (concordance) 
cosmological model only mimics the real underly ing one, which they 
call the Rh — ct Universe (see also Melia 2012 for a more popular 
description). For simplicity of notation, it will be referred to as the 
'Ht=l model' hereunder. 

As one of the motivations to introduce their model, MS 12 point 
out that the currently observed coincidence that Rh = cto (where 
to denotes the current age of the Universe) is possible within ACDM 
models only once, i.e. exactly now. On the other hand, as they ar- 
gue, such an equality [i.e. _Rh = c/H(t)] is obtained for all times 
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in a straightforward manner if one applies the Weyl postulate when 
constructing the space- time metric. 

MS 12 start by defining the gravitational radius (better known as 
Schwarzschild radius) as 

2GM(R h ) 



Rh 



(1) 



[equation (8) of MS12], where M(Rh) is the mass enclosed within 



p(t) 



(2) 



(3) 



M(R h 
and 

V - 4?F R 3 

Vprop — ^ '"'h 

is the proper volume, defined by MS 12 such that the comoving den- 
sity of particles within it remains fixed as the Universe expands. Here, 
p(t) denotes the energy density. 

MS 12 later argue that the gravitational radius is equivalent to 
the Hubble radius. We doubt that the definition of the Schwarzschild 
radius makes sense on such large scales. One of the crucial assump- 
tions in deriving the Schwarzschild radius is a static space-time. 
This condition is certainly not true when expanding regions are con- 
sidered. Nevertheless we will continue to follow the argumentation 
given in MS 12 and point out another problem in their theory. This 
is additional to what has been described in two earlier papers deal- 
ing with some flaws of the MS12 model Jvan Oirschot et all [2010; 
iLewis & van Oirschotll2012h . 

Combining the Equations Q} - © with the Friedmann equation 
for a flat Universe 
$ttG . . 



H\t) 



(4) 



one finds 



Rh{t) 



H(t) 



(5) 



which is the usual definition of the Hubble radius. Using Rh(t) = 
a(t)rh and H(t) = a(t)/a(t), we find for the comoving Hubble 
radius rh: 



rh = 



a(t) 



(6) 



The next step of MS 12 argumentation contains another flaw. They 
claim that needs to be constant because 'otherwise Vp rop would 
not represent the volume within which the particle density is constant 
in the comoving frame' . This argument is, however, not convincing. 
In a homogeneous Universe, the size of a volume - and thus the value 
of rh - does not affect the particle density within that volume. There- 
fore, there is no reason to assume that rh is constant. Consequently, 
it can neither be implied that a is constant, nor that the scale factor is 
linear in time, nor that the equation of state (EOS) w — —1/3 (see 
hereunder). 

If the argumentation of MS 12 were correct, that would have 
some odd consequences for the evolution of the Universe, which 
we will examine in the following . Being a special case of a power- 
law cosmology (see lKumadl2012l for a recent discussion, and refer- 
ences therein), in which the scale factor a(t) oc t a , this model has a 
simple mathematical form. However, it is arguably not more plausi- 
ble from the astrophysical point of view than the standard ACDM. 
Firstly, as can be seen in equation (|6}, a constant rh would im- 
ply that a is constant and thus the deceleration parameter equal to 
zero, q(z) = 0. Secondly, from Equations <[TJ — <[3j, it follows that 
Rh oc M(Rh) oc pRh, so we get that pR^ — const at all times. In 



combination with Rh(t) = a(t)rh and assuming rh to be constant, 
this leads directly to the following scaling for energy density, which 
must hold in the H t — 1 model: 



[Ht = l] p(t)oca(t)~ 



(7) 



It is easy to check that it is equivalent to require the EOS to be w = 
p/p = —1/3, as discussed in MS12. 

Recall that according to MS 12 all other homogeneous and 
isotropic cosmological models are ruled out due to (incorrect, as we 
believe) theoretical reasoning. These models would not even have 
consistent general relativistic solution. Thus, a homogeneous dust 
Universe would not have a general relativistic solution because it has 
w = 0. Neither would inflation be possible. The contents of the Uni- 
verse would need to be fine-tuned in such a way that the overall effec- 
tive EOS equals —1/3 at all times. Such a high degree of fine-tuning 
is very hard to achieve. Another problem is that the growth of fluctu- 
ations in the H t — 1 model is driven by a (negative) pressure term, 
which looks the same no matter the perturbation length. It is thus 
somewhat unclear how the large-scale structure would emerge within 
this model. 

Finally, however, neither simplicity nor plausibility are the de- 
cisive factors for a cosmological model to be accepted as a proper 
description of the physical Universe. The condition under which we 
can treat the model as valid, or at least a good approximation to re- 
ality, is that it has successfully passed confrontations with multiple 
observational data. Not all tests, applicable elsewhere, can be cur- 
rently performed within the H t = 1 model, owing to the lack of the- 
oretical framework in some domains. For example, the perturbation 
theo ry of gravitational instability, so successfully applied in ACD M 
(e.g. iDavis et alj[201ll : iNusser & Davisll201 ll : iBilicki et alj|201ll for 
recent results), could not be used in the discussed model. As shown 
in MS 12, the growth of density fluctuations is qualitatively different 
there from the standard cosmology and its quantitative details still 
await investigation. However, as we will show below, such analyses 
might already be no more than of academic interest, as the model can 
be effectively falsified with other, simpler methods. Examining the 
consequences of the H t = 1 model for the expansion history of the 
Universe and its various indicators leads us to conclude that the MS 12 
model is strongly disfavoured by observations. 



3 SUPERNOVA TYPE IA DATA 

One of the most important tests of the cosmic expansion history con- 
sists in comparing redshifts and distances of far-away objects with 
predictions of particular models. For example, in case of SNela, the 
quantity used is the luminosity distance, derived by lightcurve fitting 
from observed peak luminosities and analysed afterwards as a func- 
tion of cosmological redshift. 

Assuming spatial flatness, as was done in MS120the luminosity 
distance in ACDM models with a cosmological constant (dark energy 
equation of state w — —1) is given by 

z 

dv = -£-(l + z) [ ^ r-. (8) 

Ho K 'J [f7 m (l + C) 3 + av] 1/2 
o 

In the H t = 1 model, it takes the following closed-form expression: 
[Ht = l] d L = R h (t )(l + z)\n{l + z) , (9) 

1 The analysis can be generalized to non-flat models at the expense of more 
complicated expressions for particular quantities. 
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where Rh(to) = c/Hq. In principle, comparing this curve with avail- 
able data could be used to directly verify the validity of the model. 
Indeed, MS 12 discuss it shortly and even admit that this form of the 
luminosity dist ance does not fit t he current sample of Type la super- 
novae (see e.g. lRiess et alj|2004l) . They however argue that this dis- 
crepancy is present only for redshifts of ~ 0.5, where - as they claim 
- different SNela catalogues are inconsistent with each other. In par- 
ticular, they refer to the fact that analyses made by different authors 
do not yield a consistent value of the redshift jz acc at which, within 
ACDM, there was a transition from past deceleration to currently ob- 
served acceleration of the universal expansion. 

Putting these issues aside, we only mention a known caveat re- 
lated to using luminosity distances of SNela, or their more observer- 
friendly variant, distance moduli, as direct pr obes of cosmic expan- 
sion. When using the SALT light-curve fitter jGuv et al.ll2005|]2007t) 
some nuisance parameters that are needed to determine the distance 
moduli from the raw data, are fitted simultaneously with the cosmo- 
logical parameters. Thus, the distance muduli in the data sets are 
somewhat mo del-dependent. While the MLCS2k2 light-curve fitter 
dJha etalj [2007) avoids these issues, other systematics are present 
(e.g. the reddening parameter Ry has to be fixed by hand). Thus, al- 
though both types of light-curve fitters have some systematics, these 
are of different nature. So when accelerated expansion is found in- 
dependently by using different fitters, this should be a fairly robust 
result. 

When discussing the luminosity distance, MS 12 question the va- 
lidity of introducing into analyses such parameters as the decelera- 
tion parameter (called the acceleration parameter therein). However, 
it is by examining this parameter that one can put additional con- 
straints on the validity of the Ht=l model. It has been widely 
used for many decades as a direct probe of the cosmic expansion 
and it has straightforward interpretation. For example, within ACDM 
models (also those wit h A = 0), its cu rrent value is exactly equal to 
go = fim/2 - Ov (e.g. |Peacocklll999l) . 

In general, the deceleration parameter is defined as 



«(*) 



a a 
U 



H'( 



H(z) 



(1 + *) 



(10) 



where prime denotes differentiation with respect to redshift. The Hub- 
ble parameter is defined in the standard way, H = a/ a, where dot 
stands for time derivative. From definition then, q(z) is negative if 
the universal expansion is accelerating. In the H t = 1 model, 



[Ht = l] 







(11) 



at all times. However, various independent analyses of SNela have 
shown that the current value of the deceleration parameter is signifi- 
cantly less then zero and that it is not constant in time. This conclu- 
sion is i ndependent of the c hoice of the light-curve fitter. 

In ISeikel & Schwarzl d2009l> . a general test of accelerated ex- 
pansion has been performed using the MLCS2k2 and SALT fitters 
and different data sets. In this analysis, strong evidence for a current 
phase of accelerated expansion has been found. Another exa mple is 
the kinematic analysis performed bv lGuimaraes et all d2009h , which 
does not assume a particular underlying cosmological model. Among 
several parametrisations used therein, of particular in terest is the one 
with constant q(z). A fit to the Union Compilation (Ko walski et alj 
2008) gave q = -0.34 ± 0.05 (la errorbars). Even with the as- 
sumption of q being constant, the Ht = 1 model is at odds with 
this result. Relaxing the assumption, strongly points towards a q(z), 
which is changing with time and significantly less than zero in the 
late Universe - a result that is inconsistent with the H t = 1 model. 
This has been shown in numerous analyses. Especially those includ- 
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Figure 1. Deceleration parameter estimated from the Union2.1 SN data set. 
The dark blue line is the non-parametric smooth reconstruction based on the 
Gaussian process, with 68% and 95% confidence levels in lighter blue. The 
dashed green line presents the prediction of the H t = 1 model, whereas the 
dotted black line is for ACDM. 



ing newer observations of SNela at redshifts bigger than 1, as well 
as other probes of cosmic expansion (chronometers, baryon acoustic 
oscillations) unambiguously point to a transition from past decelera- 
tion to current acceleration (see e.g. ISantos et al.|201 ltlLu et alj201 ll ; 



iNair et alj|2012l : iGiostri et al.ll2012tlSeikel et al 11201 2al) 

As an example, we show the deceleration parameter q(z) ob- 
tained by applying Gaussian processes - a non-parametric and fully 
Bayesian approach f or data smoothing - to the distance data given 
in the Union2.1 set dSuzukietal]|2012h . This result is taken from 
Seik el et al] d2012al) , where the recently developed GaP P3 (Gaussian 
Process es in Python) package is used for the analysis. ISeikel et"al] 
(2012a) also contains an introduction to Gaussian processes. Here, 
we have added the line q = representing the Ht — 1 model for 
comparison (see Figure[T}. While the ACDM curve is well within the 
95% confidence level of the reconstructed q(z), the Ht — 1 model 
is clearly inconsistent with the data. 



4 HUBBLE RATE DATA 

We now proceed to an analysis based on observational Hubble pa- 
rameter data compiled from several sources, independent of SNela. 
We combine measurements obtained with two methods: cosmic 
chronometers (CC), which are mainly early type and/or passively 
evolving galaxies, and radial baryon acoustic oscillations (BAO) from 
galaxy clustering in redshift surveys. Using these probes should guar- 
antee robustness of our results: possible systematics are expected to 
be different in these two approaches and distinct from those in SNela 
analyses. 

Based on our dataset, shortly described later in this Section, 
we will estimate the Hubble rate and its derivatives with respect 
to redshift. For the purpose of our analysis, we will normalise the 
Hubble parameter by its current val ue: h(z) = H(z)/H q. We use 
Ho = 70.4 ± 2.5 kms^Mpc -1 jKomatsu et al]|201ll) through- 
out, although changing it to Ho = 74.2 ± 3.6 kms _1 Mpc~ 1 

2 GaPP is available for download at 

http : / / www .acgc.uct.ac.za/ ~s e ike 1/ GAPP /index . html 
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dRiess et alj|2009» , preferred by MS 12, would not influence our con- 
clusions. Subsequently, we use these estimates to determine q(z) and 
an d the so-called O m d iagnostic. These tests were recently applied 
bv ISeikel et all d2012bh against the standard model and the results 
were consistent with ACDM. Here we will use the same methods and 
(slightly extended) data set to test the Ht = 1 model. Before we 
present the data used, let us start with the theoretical framework. 

In the Ht = 1 model, because a(t) — H t and owing to the 
general relation a(z) = (1 + z) -1 (scale factor is normalised to unity 
at to), it follows that 



[Ht=l] h(z) = l + . 



h'(z) 



h"(z) = 



(12) 



Finally, the Oil' diagnostic, first introduced bv lSahni et al.| j2008). is 
defined as 



h\z) 



1 



(13) 



(1 + z) 3 

This quantity was proposed to test flat ACDM models because in that 
case it is always equal to the today's value of the matter density pa- 
rameter f2 m (to). In the H t = 1 model, Gin is variable and depends 
on redshift only: 

2 + z 



[Ht=l] Oi\\z) 



3 + 3z + z 2 



(14) 



In particular, the Ht — 1 model predicts that Oi!' (to) = 2/3. 

In order to estimate the above-d escribed quantities from obser- 
vations, we will use the same data as I Seikel et al.l J2012bh . with two 
data points added. The data set consists of 26 measurements of the 
Hubble parameter from the following sources: 

• cosmic chronom eters (CC): 18 data p o ints compiled by 
IMoresco et ail d2012bl) fro m lSimon etail d2005l) , Istern et all feoifj) 



and iMoresco et alj ( 120123) . spanning a redshift nterval of 0.09 ^ 
z < 1.75; 

• radial BAO meas ur ements: 8 da t a poin ts compiled from 
iGaztanaga etail d2009t) iBeutleret ail 1 2011 ). IChuang & Wand 
d201 ll) . iBlake et alj d2012h and iReid etail d2012h . covering redshifts 
of 0.106 < z < 0.73. 

For the reconstruction of the Hubble rate and its derivatives, we 
use again the Gaussi an process pa ckage GaPR Such an analysis has 
already been done in lSeikel et al.l d2012bl) . As we use practically the 
same sample (with only two BAO measurements added), our findings 
concerning th e observed h(z) etc . cannot be much different. How- 
ever, unlike in Seik el et all d2012bl) , here we make comparisons with 
the Ht = l model, which predicts particularly simple forms of the 
discussed functions, as given above in Equations d I lb - d!2t and d!4t 
[recall that q(z) = in the Ht=l model]. The comparisons are 
illustrated in Figures [2] and [3] presenting respectively h(z), h'(z), 
h"(z), q(z) and (z). For reference, we also plot the curves for 
the st andard flat ACDM model with fi m = 0.275 dKomatsu etail 

l2oTH) . 

The comparison with observational reconstructions is striking. 
The predictions of the Ht = 1 model are in clear disagreement with 
observations, especially at low redshifts. Very strong constraints are 
obtained here from the first derivative of the Hubble parameter. The 
Ht = 1 model requirement that h'(z) = 1 is several sigma off the 
data at redshifts z < 0.5. The same conclusion can be made from the 
reconstruction of q(z). A coasting Universe, which is predicted by the 
H t = 1 model, is inconsistent with observations. The discrepancy is 



even more pronounced in the case of the 



diagnostic (Figure 
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Figure 2. Top to bottom: normalised Hubble parameter, h(z) = H(z)/ Hq, 
and its two derivatives with respect to redshift. Top panel shows data points 
from cosmic chronometers and BAO (red crosses with errorbars). In all the 
panels, the dark blue lines are the non-parametric smooth reconstructions of 
the data, with 68% and 95% confidence levels in lighter blue; the dashed green 
lines present the predictions of the H t = 1 model, whereas the dotted black 
lines are for ACDM. 



[3}. Here, the Ht — 1 model predicts much larger values than those 
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Figure 3. The deceleration parameter q(z) (top) and the O^' diagnostic of 
cosmic acceleration (bottom). The dark blue line is the non-parametric smooth 
reconstruction based on the CC+BAO data, with 68% and 95% confidence 
levels in lighter blue. The dashed green line presents the prediction of the 
H t = 1 model, whereas the dotted black line is for ACDM. 



observed, up to z < 1.5. In our opinion, this single plot could be 
sufficient to decisively falsify the H t = 1 model. 



5 SUMMARY AND CONCLUSIONS 

• We anal ysed some of the prediction s of the cosmological model 
proposed by Meli a & Shevchukl d2012n . in which the scale factor 
growths linearly with time and the Hubble parameter is always the 
inverse of the age of the Universe (the 'Ht = 1 model'). We started 
by pointing out some inconsistencies in the derivation of the model 
and its bizarre consequences for the nature of the Universe. 

• We referred to several analyses of supernova la data that are 
all inconsistent with the Ht — 1 model. We presented a plot, 
which compares the deceleration parameter predicted by this model 
[q(z) = 0] with a model-independent reconstruction of q(z) obtained 
by applying Gaussian processes to the Union2.1 set. 

• We used Hubble parameter data obtained from cosmic 
chronometers ( 18 da ta points with 0.09 ^ z 1.75), gathered by 
Mores co et al. I d2012bl) . and from baryon acoustic oscillations (8 mea- 
surements in total), combined from several sources, that spanned a 
redshift interval 0.106 ^ z ^ 0.73. We compared the observational 



data with the predictions of the H t = 1 model regarding the fol- 
lowing parameters: the Hubble parameter H(z), its derivatives with 
respect to redshift H'(z) and H"[z), as well as the O m (z) diagnostic 
and the deceleration parameter q(z). In order to obtain derivatives of 
functions base d on discrete data, w e applied the recently developed 
GaPP package dSeikel et alj|2012<j) . 

• For all these diagnostics we find severe discrepancies between 
the Ht = 1 model and observations, especially at low redshifts. 
The most pronounced disagreements are for the first derivative of the 
Hubble parameter, the O m diagnostic and the deceleration parame- 
ter. The analysed model completely fails to reproduce their observed 
behaviou r for z < 0.5. This strong ly suggests that the cosmological 
model of lMelia & Shevchukl J2012h is not a proper description of our 
Universe. 

• The Rh = ct model proposed by iMelia & Shevchukl d2012l) is 
mathematically simple, but it does not describe the Universe we live 
in. 
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